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FIELD OF THE INVENTION 
This invention relates generally to colonoscopy techniques and more particularly to a 
system for processing the image of a bowel to remove bowel contents from the image and for 
detecting polyps in the digital image. 

BACKGROUND OF THE INVENTION 
As is known in the art, a colonoscopy refers to a medical procedure for examining a colon 
to detect abnormalities such as polyps, tumors or inflammatory processes in the anatomy of the 
colon. The colonoscopy is a procedure which consists of a direct endoscopic examination of the 
colon with a flexible tubular structure known as a colonoscope which has fiber optic or video 
recording capabilities at one end thereof The colonoscope is inserted through the patient's anus 
and directed along the length of the colon, thereby permitting direct endoscopic visualization of 
colon polyps and tumors and in some cases, providing a capability for endoscopic biopsy and 
polyp removal. Although colonoscopy provides a precise means of colon examination, it is time- 
consuming, expensive to perform, and requires great care and skill by the examiner. And, the 
procedure also requires thorough patient preparation including ingestion of purgatives and 
enemas, and usually a moderate anesthesia. Moreover, since colonoscopy is an invasive 
procedure, there is a significant risk of injury to the colon and the possibility of colon perforation 
and peritonitis, which can be fatal. 
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To overcome these drawbacks, the virtual colonoscopy was conceived. A virtual 
colonoscopy makes use of images generated by computed tomography (CT) imaging systems 
(also referred to as computer assisted tomography (CAT) imaging systems). In a CT (or CAT) 
imaging system, a computer is used to produce an image of cross-sections of regions of the 
5 human body by using measure attenuation of X-rays through a cross-section of the body. In a 
virtual colonoscopy, the CT imaging system generates two-dimensional images of the inside of an 
intestine. A series of such two-dimensional images can be combined to provide a three- 
dimensional image of the colon. While this approach does not require insertion of an endoscope 
into a patient and thus avoids the risk of injury to the colon and the possibility of colon 
1 0 perforation and peritonitis, the approach still requires thorough patient preparation including 

purgatives and enemas. Generally, the patient must stop eating and purge the bowel by ingesting 
(typically by drinking) a relatively large amount of a purgative. Another problem with the virtual 
*g: colonoscopy approach is that, the accuracy of examinations and diagnosis using virtual 
SI colonoscopy techniques is not as accurate as is desired. This is due, at least in part, to the 
Ife relatively large number of images the examiner (e.g. a doctor) must examine to determine if a 
OR P°lyp> tumor or an abnormality exists in the colon. 

y It would, therefore, be desirable to provide a virtual colonoscopy technique which 

p removes the need for bowel cleansing. It would also be desirable to provide a virtual colonoscopy 
0 technique which removes the need for thorough patient preparation. 

It would further be desirable to provide a technique which increases the accuracy of 
examinations and diagnosis using virtual colonoscopy. It would be further desirable to provide a 
technique which reduces the number of images an examiner (e.g. a doctor) must examine to 

2 5 determine if a polyp, tumor or an abnormality exists in the colon. It would be further desirable to 

provide a technique for automatic detection of polyps, tumors or other abnormalities in the colon. 

SUMMARY OF THE INVENTION 
In accordance with the present invention, an apparatus for performing a virtual 

3 0 colonoscopy includes a system for generating digital images, a storage device for storing the 
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digital images, a digital bowel subtraction processor coupled to receive images of a colon from 
the storage device and for processing the received digital images of the colon to digitally remove 
the contents of the colon from the image. With this particular arrangement, a system which 
provides accurate results without the need for thorough patient preparation (e.g. without the need 
5 for bowel cleansing) is provided. The digital bowel subtraction processor (DBSP) receives image 
data from the image database and processes the image data to digitally remove the contents of the 
bowel from the digital image. The DBSP can then store the image back into the image database. 
Since the DBSP digitally subtracts the contents of the bowel, the patient undergoing the virtual 
colonoscopy need not purge the bowel in the conventional manner which is know to be 

1 0 unpleasant to the patient. The system can further include an automated polyp detection processor 
coupled to receive images of a colon from the storage device and for processing the received 

^ digital images of the colon to detect polyps in the colon image. 

fli The automated polyp detection processor (APDP) receives image data from the image 

3S storage device (which may be provided as an image database, for example) and processes the 
ffl image data to detect and /or identify polyps, tumors, inflammatory processes, or other 
E ~ irregularities in the anatomy of the colon. The APDP can thus pre-screen each image in the 
% database such that an examiner (e.g. a doctor) need not examine every image but rather can focus 
O attention on a subset of the images possibly having polyps or other irregularities. Since the CT 
£§ 0 system generates a relatively large number of images for each patient undergoing the virtual 
s ~* colonoscopy, the examiner is allowed more time to focus on those images in which it is most 
likely to detect a polyp or other irregularity in the colon. The APDP can process images which 
have been generated using either conventional virtual colonoscopy techniques (e.g. techniques in 
which the patient purges the bowel prior to the CT scan) or the APDP can process images in 

2 5 which the bowel contents have been digitally subtracted (e.g. images which have been generated 

by DBSP). 

In accordance with a still further aspect of the present invention a technique for digital 
bowel subtraction includes the steps of applying a threshold function to the image data and 

3 0 selecting all image elements above the threshold for further processing, applying a gradient to the 
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original image to identify "shoulder regions" and then performing an image dilation step on the 
shoulder region (i.e. an expansion is performed on selected pixels in the shoulder region). The 
bowel contents are digitally subtracted based on the gradient analysis. Thus, after selecting and 
dilating the gradient the subtraction of the bowel contents is performed. Mathematically this can 
5 be represented as Subtracted Image Values ~ Original Image Values- (Threshold Values + 
Gradient Values). 

In accordance with a still further aspect of the present invention a technique for digital 
bowel subtraction includes the steps of scanning across a matrix of digital values which represents 
10 an image, identifying regions corresponding to regions of a colon, identifying regions 

corresponding to regions of air and bowel contents in the colon and subtracting one region from 
^ the other to provide an image with the bowel contents removed. With this technique, a raster 
ypt based searching method is provided. By scanning an image in a raster pattern and applying 
ry threshold values in a predetermined logic sequence, pixels representing air which are located 
i J 5 proximate pixels representing bowel wall are found. Once a region corresponding to a boundary 
111 region between bowel contents and a bowel wall is found, the pixels which represent the bowel 
1~ contents can be subtracted from the image. In one embodiment, the pixels representing the bowel 

contents are subtracted from the image by setting the values of the pixels to value corresponding 
O to air. 

3° 

p= In accordance with a still further aspect of the present invention a technique for automatic 

polyp detection includes the steps of generating a polyp template by obtaining a CT image which 
includes a polyp and excising the polyp from the image and using the excised polyp image as the 
template. Next, a portion of a CT image is selected and a polyp identification function is applied 
25 to the selected portion of the CT image. Next, a correlation is performed between the output of 
the identification function and the image. With this particular arrangement, a technique for 
automated polyp detection is provided. 

In accordance with a still further aspect of the present invention, a second technique for 
3 0 automatic polyp detection includes the steps of moving a test element along a boundary of a 
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bowel When the test element is rolled around the bowel perimeter, the changes in direction of 
the test element are detected as the changes occur. The features of the bowel perimeter are 
classified based upon the path of the test element. Thus, an advance knowledge of the geometric 
characteristic of the polyp or other irregularity being detected is required. The process of 
5 classifying the features of the bowel are accomplished by marking three points where turns exist. 
The points are selected by looking at the changes in slope (e.g. the derivative). Next the distances 
between the marked points are computed and the ratios between the distances are used to identify 
bowel features or characteristics of the bowel In one embodiment, the test element is provided 
having a circular shape and thus is referred to as a virtual rolling ball, a virtual ball, or more simply 

10 a ball- 

In accordance with a still further aspect of the present invention, a third technique for 
J3 automatic detection of structures, including but not limited to polyps, includes the steps of 

11 applying a template to a segmented bowel image corresponding to the bowel perimeter and 

835 computing the distances between points on the template perimeter and the perimeter bowel points 

01 present within a window. Next, a determination is made as to whether the distances are equal 
When the template becomes centered within a lesion, then the distances from one or more points 

5 on the template to bowel boundary points of the structure (e.g. a lesion)become substantially 
O equal. The distances may be measured from a number of points on the template. For example, a 
;2jj0 center point of the template, perimeter points of the template or other points on the template may 
N 1 be used. The particular points to use on the template are selected in accordance with a variety of 
factors including but not limited to the template shape and the physical characteristics (e.g. shape) 
of structure being detected. To make a determination of when the computed distances are 
substantially equal, a standard deviation of a group of those distances can be computed. The 

2 5 template location at which the standard deviation values approach a minimum value corresponds 

to a location at which a structure having a shape similar to the template shape exists. 
Alternatively, an average distance from a point on the template (e.g. a template center point) to 
the perimeter of the boundary structure (e.g. a bowel wall) can be computed. In this case, a 
location at which a structure having a shape similar to the template shape exists can be found 
30 when the average distance from the point on the template to the boundary structure perimeter 



points reaches a minimum value. The above process can be carried out in one plane or in three 
orthogonal planes. The point where the standard deviation of a group of those distances 
approaches a minimum value is the point which should be marked as a center of a suspected 
lesion. In the case where the technique is run in three orthogonal planes, those lesions that were 
tagged in 2 of 3 or 3 of 3 planes can be finally tagged as suspicious regions. This technique can 
thus be used to distinguish a fold from a polyp. With this particular arrangement, a technique 
which searches for patterns of distance is provided. The technique is thus relatively 
computationally expensive but it is also relatively rigorous. One advantage of this technique is 
that does not matter how large the template is relative to the polyp so one scanning works for all 
lesions. 

BRIEF DESCRIPTION OF THE DRAWINGS 
The foregoing features of the invention, as well as the invention itself may be more fully 
understood from the following detailed description of the drawings, in which: 

FIG. 1 is a block diagram of a system for digital bowel subtraction and automatic polyp 
detection; 

FIGs. 1 A - ID a series of views illustrating a digital image of a bowel before and after 
processing via a digital bowel subtraction processor; 

FIG. 2 is a flow diagram showing the steps in a virtual colonoscopy performed using 
digital bowel subtraction; 

FIG. 3 is a flow diagram showing the steps in a first method for performing digital bowel 
subtraction; 

FIG. 4 is a diagram showing the volume averaging area in a portion of a colon; 
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FIGs. 4 A - 4C are plots of pixel values vs. pixel locations in the formation of a bowel 
wall/air boundary; 

FIG. 5 is a flow diagram showing the steps in a second method for performing digital 
5 bowel subtraction; 

FIGs. 5 A is diagrammatic view of an image of the type generated by a CT system and 
having a local window disposed thereover; 

10 FIG. 5B is a portion of the image shown in FIG. 5 A taken along lines 5B-5B showing a 

polyp obscured by bowel contents; 

yg FIG. 5C is a portion of the image shown in FIG. 5 A taken along lines 5B-5B after the 

jJt digital bowel subtraction process and in which the polyp is clearly visible; 

|5 

m Fig. 5 D is a pixel analysis map; 

y FIG. 6 is a flow diagram showing the steps to digital subtract bowel contents in a region 

O between opacified and non-opacified material in a bowel; 

^ FIGs. 6A and 6B are a series of plots illustrating the threshold process; 

FIG 7 is a flow diagram showing the steps of a template matching technique to 
automatically detect a polyp in a CT image of a bowel; 

25 

FIG. 8 is a flow diagram showing the steps of a rolling ball technique to automatically 
detect a polyp in a CT image of a bowel; and 



FIGs. 8A-8H are examples of automatic polyp detection in a CT image of a bowel using 
3 0 the rolling ball technique; 
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FIG. 81 is an example of automatic polyp detection using the rolling ball technique in a 
three dimensional image; 

5 FIG. 9 is a flow diagram showing the steps of a distance matching technique to 

automatically detect a structure in a CT image of a bowel; and 

FIGs. 9A-9G are examples of automatic polyp detection in a CT image of a bowel using 
the distance matching technique. 

10 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
Before describing a virtual colonoscopy system which includes a digital bowel subtraction 
m processor (DBSP) and / or automated polyp detection processor (APDP) and the operations 
;SJ performed to digital cleanse a bowel and automatically detect a polyp, some introductory 
ffe concepts and terminology are explained. 

^ A computed tomography (CT) system generates signals which can be stored as a matrix of 

digital values in a storage device of a computer or other digital processing device. As described 
Q herein, the CT image is divided into a two-dimensional array of pixels, each represented by a 
So digital word. One of ordinary skill in the art will recognize that the techniques described herein 
^ are applicable to various sizes and shapes of arrays. The two-dimensional array of pixels can be 
combined to form a three-dimensional array of pixels. The value of each digital word corresponds 
to the intensity of the image at that pixel. Techniques for displaying images represented in such a 
fashion, as well as techniques for passing such images from one processor to another, are known. 

25 

As also described herein, the matrix of digital data values are generally referred to as a 
"digital image" or more simply an "image" and may be stored in a digital data storage device, such 
as a memory for example, as an array of numbers representing the spatial distribution of energy at 
different wavelengths in a scene. 

30 
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Each of the numbers in the array correspond to a digital word typically referred to as a 
"picture element" or a "pixel" or as "image data." The image may be divided into a two 
dimensional array of pixels with each of the pixels represented by a digital word. Thus, a pixel 
represents a single sample which is located at specific spatial coordinates in the image. 

5 

It should be appreciated that the digital word is comprised of a certain number of bits and 
that the techniques of the present invention can be used on digital words having any number of 
bits. For example, the digital word may be provided as an eight-bit binary value, a twelve bit 
binary value, a sixteen but binary value, a thirty-two bit binary value, a sixty-four bit binary value 
10 or as a binary value having any other number of bits. 

It should also be noted that the techniques described herein may be applied equally well to 
yn either grey scale images or color images. In the case of a gray scale image, the value of each 
St digital word corresponds to the intensity of the pixel and thus the image at that particular pixel 
tffe location. In the case of a color image, reference is sometimes made herein to each pixel being 
fjft represented by a predetermined number of bits (e.g. eight bits) which represent the color red (R 

bits), a predetermined number of bits (e.g. eight bits) which represent the color green (G bits) and 
Q a predetermined number of bits (e.g. eight bits) which represent the color blue (B-bits) using the 
Q so-called RGB color scheme in which a color and luminance value for each pixel can be computed 
2[0 from the RGB values. Thus, in an eight bit color RGB representation, a pixel may be represented 
N= by a twenty-four bit digital word. 

It is of course possible to use greater or fewer than eight bits for each of the RGB values. 
It is also possible to represent color pixels using other color schemes such as a hue, saturation, 
25 brightness (HSB) scheme or a cyan, magenta, yellow, black (CMYK) scheme. It should thus be 
noted that the techniques described herein are applicable to a plurality of color schemes including 
but not limited to the above mentioned RGB, HSB, CMYK schemes as well as the Luminosity 
and color axes a & b (Lab) YUV color difference color coordinate system, the Karhunen-Loeve 
color coordinate system, the retinal cone color coordinate system and the X, Y, Z scheme. 

30 
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Reference is also sometimes made herein to an image as a two-dimensional pixel array. 
An example of an array size is size 512 x 512, One of ordinary skill in the art will of course 
recognize that the techniques described herein are applicable to various sizes and shapes of pixel 
arrays including irregularly shaped pixel arrays. 

An "image region" or more simply a "region" is a portion of an image. For example, if an 
image is provided as a 32 X 32 pixel array, a region may correspond to a 4 X 4 portion of the 32 
X 32 pixel array. 

In many instances, groups of pixels in an image are selected for simultaneous 
consideration. One such selection technique is called a "map" or a "local window," For example, 
if a 3 X 3 subarray of pixels is to be considered, that group is said to be in a 3 X 3 local window. 
One of ordinary skill in the art will of course recognize that the techniques described herein are 
applicable to various sizes and shapes of local windows including irregularly shaped windows. 

It is often necessary to process every such group of pixels which can be formed from an 
image. In those instances, the local window is thought of as "sliding" across the image because the 
local window is placed above one pixel, then moves and is placed above another pixel, and then 
another, and so on. Sometime the "sliding" is made in a raster pattern. It should be noted, 
though, that other patterns can also be used. 

It should also be appreciated that although the detection techniques described herein are 
described in the context of detecting polyps in a colon, those of ordinary skill in the art should 
appreciate that the detection techniques can also be used search for and detect structures other 
than polyps and that the techniques may find application in regions of the body other than the 
bowel or colon. 

Referring now to FIG. 1, a system for performing virtual colonoscopy 10 includes a 
computed tomography (CT) imaging system 12 having a database 14 coupled thereto. As is 
known, the CT system 10 produces two-dimensional images of cross-sections of regions of the 
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human body by measuring attenuation of X-rays through a cross-section of the body. The images 
are stored as digital images in the image database 14, A series of such two-dimensional images 
can be combined using known techniques to provide a three-dimensional image of the colon. A 
user interface 16 allows a user to operate the CT system and also allows the user to access and 
5 view the images stored in the image database. 

A digital bowel subtraction processor (DBSP) 18 is coupled to the image database 14 and 
the user interface 16. The DBSP receives image data from the image database and processes the 
image data to digitally remove the contents of the bowel from the digital image. The DBSP can 

1 0 then store the image back into the image database 14. The particular manner in which the DBSP 
processes the images to subtract or remove the bowel contents from the image will be described 
in detail below in conjunction with Figs. 2-6. Suffice it here to say that since the DBSP digitally 

y3 subtracts or otherwise removes the contents of the bowel from the image provided to the DBSP, 

J the patient undergoing the virtual colonoscopy need not purge the bowel in the conventional 

§6 manner which is know to be unpleasant to the patient. 

The DBSP 1 8 may operate in one of at least two modes. The first mode is referred to a 
y raster mode in which the DBSP utilizes a map or window which is moved in a predetermined 
O pattern across an image. In a preferred embodiment, the pattern corresponds to a raster pattern. 
? 3[0 The window scans the entire image while threshold values are applied to pixels within the image 
^ in a predetermined logic sequence. The threshold process assesses whether absolute threshold 

values have been crossed and the rate at which they have been crossed. The raster scan approach 
looks primarily for "air" pixels proximate (including adjacent to) bowel pixels. The processor 
examines each of the pixels to locate native un-enhanced soft tissue, As a boundary between soft 
25 tissue (e.g. bowel wall) and bowel contents is established, pixels are reset to predetermined values 
depending upon which side of the boundary on which they appear. 
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The second mode of operation for the DBSP 18 is the so-called gradient processor mode. 
In the gradient processor mode, a soft tissue threshold (ST) value, an air threshold (AT) value and 
a bowel threshold (BT) value are selected. A first mask is applied to the image and all pixels 
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having values greater than the bowel threshold value are marked. Next, a gradient is applied to 
the pixels in the images to identify pixels in the image which should have air values and bowel 
values. The gradient function identifies regions having rapidly changing pixel values. From 
experience, one can select bowel/air and soft tissue/air transition regions in an image by 
5 appropriate selection of the gradient threshold. The gradient process uses a second mask to 

capture a first shoulder region in a transition region after each of the pixels having values greater 
than the BT value have been marked. 

Once the DBSP 18 removes the bowel contents from the image, there exists a relatively 
1 0 sharp boundary and gradient when moving from the edge of the bowel wall to the "air" of the 
bowel lumen. This is because the subtraction process results in all of the subtracted bowel 
contents having the same air pixel values. Thus, after the subtraction, there is a sharp boundary 
yS and gradient when moving from the edge of the bowel wall to the "air" of the bowel lumen. In 
S this context, "air" refers to the value of the image pixels which been reset to a value 
|£ corresponding to air density. If left as is, this sharp boundary (and gradient) end up inhibiting the 
f§i 3D endoluminal evaluation of the colon model since sharp edges appear as bright reflectors in the 
^ model and thus are visually distracting. 

Q A mucosa insertion processor 19 is used to further process the sharp boundary to lesson 

2jo the impact of or remove the visually distracting regions. The sharp edges are located by applying a 
gradient operator to the image from which the bowel contents have been extracted. The gradient 
operator may be similar to the gradient operator used to find the boundary regions in the gradient 
subtractor approach described herein. The gradient threshold used in this case, however, typically 
differs from that used to establish a boundary between bowel contents and a bowel wall. 

25 

The particular gradient threshold to use can be empirically determined. Such empirical 
selection may be accomplished, for example, by visually inspecting the results of gradient selection 
on a set of images detected under similar scanning and bowel preparation techniques and adjusting 
gradient thresholds manually to obtain the appropriate gradient (tissue transition selector) result. 

30 
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The sharp edges end up having the highest gradients in the subtracted image. A 
constrained gaussian filter is then applied to these boundary (edge) pixels in order to "smooth" the 
edge. The constraint is that the smoothing is allowed to take place only over a predetermined 
width along the boundary. The predetermined with should be selected such that the smoothing 
5 process does not obscure any polyp of other bowel structures of possible interest. In one 
embodiment the predetermined width corresponds to a width of less than ten pixels. In a 
preferred embodiment, the predetermined width corresponds to a width in the range of two to five 
pixels and in a most preferred embodiment, the width corresponds to a width of three pixels. The 
result looks substantially similar and in some cases indistinguishable from the natural mucosa seen 
10 in untouched bowel wall, and permits an endoluminal evaluation of the subtracted images. 

_ Also coupled between the image database 14 and the user interface 16 is an automated 

Cf P°lyp detection processor (APDP) 20. The APDP 20 receives image data from the image 

database and processes the image data to detect and /or identify polyps, tumors, inflammatory 
305 processes, or other irregularities in the anatomy of the colon. The APDP 20 can thus pre-screen 
m each image in the database 14 such that an examiner (e.g. a doctor) need not examine every image 
J* but rather can focus attention on a subset of the images possibly having polyps or other 
O irregularities. Since the CT system 10 generates a relatively large number of images for each 
-q patient undergoing the virtual colonoscopy, the examiner is allowed more time to focus on those 

images in which the examiner is most likely to detect a polyp or other irregularity in the colon. 
H The particular manner in which the APDP 20 processes the images to detect and /or identify 

polyps in the images will be described in detail below in conjunction with Figs. 7-9. Suffice it here 
to say that the APDP 20 can be used to process two-dimensional or three-dimensional images of 
the colon. It should also be noted that APDP 20 can process images which have been generated 
25 using either conventional virtual colonoscopy techniques (e.g. techniques in which the patient 
purges the bowel prior to the CT scan) or the APDP 20 can process images in which the bowel 
contents have been digitally subtracted (e.g. images which have been generated by DBSP 18). 

It should also be appreciated that polyp detection system 20 can provide results generated 
3 0 thereby to an indicator system which can be used to annotate (e.g. by addition of a marker, icon 
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or other means) or otherwise identify regions of interest in an image (e.g. by drawing a line 
around the region in the image, or changing the color of the region in the image) which has been 
processed by the detection system 20. 

Referring now to FIGs. 1 A - ID in which like elements are provided having like reference 
designations throughout the several views, a series of images 19a-19d illustrating a bowel before 
(Figs. 1A, IB) and after (Figs. 1C, ID) processing via the DBSP 18 (Fig. 1) are shown. 

In Fig. 1A, the image 19a corresponds to a single slice of a CT scan. The image 19a 
includes a plurality of opacified bowel contents 20 (shown as light colored regions) and regions of 
air 22 (shown as dark colored regions). The image 19a also includes regions 24 which 
correspond to native soft tissue and regions 26 which correspond to portions of a bowel wall. 

Image 19b (Fig. IB) is three-dimensional view of the colon formed from a series of single 
slice CT scans. The opacified bowel contents 20 are shown. The opacified bowel contents 20 
impairing the view of the bowel wall and thus limit the ability to detect regions of interest along 
the bowel wall. 

In Fig. 1C, single slice CT scan image 19c is the same as image 19a except that the bowel 
contents present in the image 19a (Fig. 19 A) have been digitally subtracted from the image 19a. 
The image 19c does not include the opacified bowel contents present in Fig. 19 A. Rather the 
bowel contents 20 (Fig. 19 A) have been replaced by regions 22 (shown as dark colored regions) 
corresponding to air. The image 19c also includes the native soft tissue and bowel wall regions 
24, 26 respectively. 

As discussed above, the DBSP 18 receives an image of a section of a bowel and digitally 
subtracts the contents 20 of the bowel section. Here the operation is performed on a two- 
dimensional image and a series of such two-dimensional images (each of which has been 
processed by the DBSP 18) can be combined to provide a three-dimensional image of the bowel 
section. Such two-dimensional and three-dimensional images can then be processed via the 



15 

APDP 20 (Fig. 1) to detect polyps, tumors, inflammatory processes, or other irregularities in the 
anatomy of the colon. 

FIGs, 2-9 are a series of flow diagrams showing the processing performed by a processing 
5 apparatus which may, for example, be provided as part of a virtual colonoscopy system 10 such as 
that described above in conjunction with FIG. 1 to perform digital bowel subtraction and 
automated polyp detection. The rectangular elements (typified by element 30 in FIG. 2), herein 
denoted "processing blocks/' represent computer software instructions or groups of instructions. 
The diamond shaped elements (typified by element 64 in FIG 5), herein denoted "decision 
1 0 blocks," represent computer software instructions, or groups of instructions which affect the 
execution of the computer software instructions represented by the processing blocks. 

J3 Alternatively, the processing and decision blocks represent steps performed by functionally 

DP- 

m equivalent circuits such as a digital signal processor circuit or an application specific integrated 
|fc circuit (ASIC). The flow diagrams do not depict the syntax of any particular programming 
ffl language. Rather, the flow diagrams illustrate the functional information one of ordinary skill in 
r* the art requires to fabricate circuits or to generate computer software to perform the processing 

required of the particular apparatus. It should be noted that many routine program elements, such 
Q as initialization of loops and variables and the use of temporary variables are not shown. It will be 
iSO appreciated by those of ordinary skill in the art that unless otherwise indicated herein, the 
^ particular sequence of steps described is illustrative only and can be varied without departing from 

the spirit of the invention. 

Turning now to FIG. 2, the steps in a virtual colonoscopy are shown. As shown in step 
25 30, the virtual colonoscopy process begins by placing a contract agent in the region of the colon 
in which the CT scan will be performed. Typically, the patient ingests the contract agent. It 
should be appreciated however, that any technique for placing the contrast agent in the bowel 
may also be used. The contrast agent may be taken in small amounts with meals beginning 
approximately 48 hours or so prior to a scheduled CT exam. The contrast agent can be of any of 
30 the commercially available types such as GastrografFm, Barium or Oxilan, for example. 
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Next, as show in step 32, the CT exam takes place and images of a body region (e.g. an 
entire abdomen) are generated as shown in step 34. 

5 The image data is then sent to an analysis system and analyzed as shown in step 36. One 

or more images are selected for analysis (can be a single slice or a series of slices). This can be 
accomplished using commercially available systems such as the Vitrea System available through 
Vital Images, Inc. Minneapolis, MN. 

1 0 Then, as shown in step 38, the contents of the bowel are digitally subtracted from the CT 

images. The digital subtraction step can be performed using either of the techniques described 
below in conjunction with FIGs. 3 and 5 below. 

?! After the digital subtraction step, the results of the bowel images having the contents 

IB5 thereof digitally removed are displayed as shown in step 40. 

Referring now to FIG. 3, a first technique for digitally subtracting the contents of a bowel 
□ utilizes functions from the so-called "toolbox" found in the MATLAB computer program 
r4 (available through The Math Works, Natick, Massachusetts). This technique begins by applying a 
%Q threshold function to the image data and selecting all image elements above the threshold. 

Next as shown in step 44, a gradient is applied to the original image to identify "shoulder 
regions." One the should regions are identified, an image dilation step is performed on the 
shoulder region (i.e. an expansion is performed on selected pixels in the shoulder region). The 
2 5 gradient functions can be provided as any number of functions including but not limited to the 

Roberts or Sobel gradient functions. Generally, in those regions in which the air water interfaces 
are horizontal, it is desirable to use a horizontal gradient. 
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Next, as shown in step 46, the bowel contents are digitally subtracted based on the 
gradient analysis. Thus, after selecting and dilating the gradient, the subtraction of the bowel 
contents is performed. Mathematically this can be represented as: 

Subtracted Image Values ~ Original Image Values- (Threshold Values + Gradient Values) 

in which : 

Original Image Values = the values of the digital image as measured by the CT System and 

massaged as necessary but prior to DBS processing 

Threshold Values = a selected threshold value 

Gradient Values = values resultant from the selected gradient function 

Referring briefly to FIGs. 4 - 4C, an image of a bowel portion 49 (Fig. 4) includes an 
unopacified region 50 (Fig. 4), and an opacified region 53 (Fig. 4) and a volume averaging region 
comprised of sections 51, 52 (Fig. 4). The contrast agent ingested increases the contrast between 
the bowel contents (represented by opacified region 53) and air (represented by the unopacified 
region 50). The regions 51, 52 correspond to a boundary area between the unopacified region 53 
and the opacified region 53 . That is, a clearly defined boundary does not exist between the 
unopacified and opacified regions 49, 53. 

Referring now to FIG. 4A ? a plot 55 of pixel values in Hounsfield Units (HU) vs. pixel 
locations in a CT image reveals a first region 55a corresponding to a soft tissue region, a second 
region 55b corresponding to an opacified bowel contents region and a third region 55c 
corresponding to an air region. A transition region thus exists between the soft tissue region 55a 
and the air region 55c. . A gradient function applied to curve 55 produces boundary 57 shown in 
Fig. 4B. Thus, as shown in Fig. 4B, after application of a gradient function, and subtraction, a 
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relatively sharp transition 57a exists between regions 55a and 55c. It is this transition on which 
the mucosal insertion processor 19 operates as described above in conjunction with Fig. 1. 

As shown in Fig. 4C ? after application of the mucosal insertion process, a transition 57a 
between the soft tissue region 55a and the air region 55c is provided. 

Referring now to FIG. 5, a second technique for digitally subtracting the contents of a 
bowel includes the step of generating a pixel analysis map (PAM) from the image as show in step 
60. In one particular embodiment, the image is provided as a 5 12 X 5 12 image and the PAM is 
provided as a 7 X 8 matrix generated from the image. 

The PAM is moved across the image in a raster pattern. As the PAM is moved from 
location to location across the different regions of the image, steps 62-79 are performed. As 
shown in step 62, the central elements (upper left (ul), upper right (ur) lower left (11), lower right 
(lr) shown in Fig. 5 A) of the PAM are the focused of the analysis. Steps 64 and 66 implement a 
loop in which a search is performed until an air region is found. 

Once the air region is found, processing proceeds to step 68 where the area around the air 
element is searched to locate a bowel region as shown in step 70. If the region corresponds to a 
bowel region, then processing proceeds to step 74 where it is determined if the bowel region is 
adjacent to an air region. If the bowel region is adjacent to the air region, then the bowel contents 
are subtracted as shown in step 76. 

If in step 70 a decision is made that the selected PAM does not correspond to bowel, then 
steps 77 and 78 implement a loop in which a new PAM is selected until a bowel region is found. 

This process is repeated for each of the ul, ur, 11, lr in the local window during the raster 
pattern scan. 
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Referring to Figs. 5 A-5C in which like element are provided having like reference 
designations throughout the several views, a 512 X 512 image 80 of the type generated by a CT 
system for example, is shown having a 7 x 8 local window 82 disposed thereover. The local 
window 82 is moved across the image in a raster pattern designated as reference numeral 84. The 
5 image 80 is that of a two-dimensional section of a bowel 86. Portions of the bowel 86 have 

contents 88a (indicated by cross-hatching) and portions 88b correspond to air. The bowel 86 has 
a fold 89 therein and a plurality of polyps 90a - 90b. 

As may be more clearly seen in Fig. 5B, in which like elements in Fig. 5 A are shown 
1 0 having like reference designations, polyp 90b is obscured by the bowel contents 88a and thus may 
be difficult to see in an image generated by a CT scan. In particular, such a polyp may be 
particularly difficult to see in a three-dimensional (3D) view. After the digital subtraction process 
^ of the present invention, the image shown in FIG. 5C results. In FIG. 5C, the bowel contents 
^ have been removed and the polyp 90b is clearly visible. Thus, when a visual examination of the 
£05 CT image is conducted, the polyp 90b is exposed and can be easily viewed. 

Referring now to Fig. 5D, a pixel analysis map (PAM) 82' includes four central elements 
;3 90a-90. As described above, axial digital bowel cleansing can be performed using a marching 
5 squares routine to scan through selected regions of an image (e.g. image 80 in Fig. 5 A) to reset 
JgO pixels that fall above a bowel threshold value. The bowel threshold value may be set by the user, 
M One technique uses an air threshold and a soft tissue threshold (both of which are expressed in 

Hounsfield units) to also segment the bowel mucosal boundary, for use in polyp detection. 

Segmented pixels (the bowel boundary) are set to a value of 2500. 

25 In the technique, the central four elements 91a-91d of the PAM 82' are used as the 

marching square. An examination of the image to which the PAM 82' is applied is made by 
looking out to the far boundaries created by the "L" and "R" level pixels. It should be noted that 
in Fig. 5D> an "o" designates boundary pixels while an "x" designates outer boundary pixels. 
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In the technique one of the central pixels is selected (e.g. ul 91a) and one direction is 
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selected (e.g. the SE direction designated by reference line 92. Next, the furthest point away in 
that direction is examined. In this example, the furthest point corresponds to the point designated 
SE is selected. 

5 The difference between the value of the ul pixel 9 1 a and the value of the pixel at the 

location SE is computed. If the difference between the values of the SE and ul is large enough, 
then one can conclude that a substantial threshold exists. Then points which are closer to the ul 
pixel than the SE pixel are examined (e.g. pixels at the IL and LI locations). This process is then 
repeated for each direction. Once all the directions have been tested, then a new central pixel is 

0 selected (e.g. one of pixels 91b-91d) and the process is repeated. 

Referring now to FIG. 6, the process of searching around the pixel ul in the local window 
begins with steps 100 and 102 in which a PAM and one of the elements ul, ur, 11, lr are selected. 
In the example discussed in FIG. 6, the pixel in location ul is selected as shown in step 102. Thus 
the processing starts at pixel ul and each of the neighbors to ul (e.g. ur, 11, lr) are examined. It 
should be appreciated, however, that the process could also start one of the other pixels (i.e. ur, 11, 
lr) and then the corresponding neighbors would be examined. 

Considering first neighbor ur, decision is made as to whether ur is above the bowel 
?0 threshold as shown in step 106. This process is explained further in conjunction with FIG. 6 A. If 
the value of ur is above the bowel threshold, then processing flows to step 108 in which the pixel 
value is subtracted to reset the value to the air value. Processing then proceed to step 114 m 
which the next neighboring pixel is examined. 

2 5 If in step 106 it is determined that the value of ur is not above the bowel threshold, then a 

determination is made as to if any more pixel exist to the right of ur as shown in step 110. If such 
pixels do exist then the system looks one pixel further to the right and processing returns to step 
106. If no such pixels exist, then processing proceeds to again to step 1 14. 
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The above steps are repeated for each of the pixels ul, 11, lr as indicated in step 116. 
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Referring now to FIGs. 6A and 6B, the threshold process is illustrated. There are two 
parts to the threshold process. The first part is to determine the air-bowel boundary. The second 
part requires a decision as to how to treat shoulder region 120. 

5 

If the difference between the extreme neighbor uE and the close neighbor PI is above a 
gradient threshold, then use that as well as an indication of the air-bowel interface and all values 
between the starting pixel (i.e. ul in the present example) and the far extreme pixel (i.e. uE in the 
present example and which is the bowel boundary) are reset to air values. That is, the values for 

1 0 pixels P 1 , P2 are reset to provide the boundary 1 22 shown in FIG. 6B . The resultant pixel values 
(i.e. the reset values for pixels PI, P2) are the values upon which the subtraction is based. The 
pixel values are preferably set in this way since it is computationally efficient to trigger with the 

Sj smallest number and because the smallest polyps of interest are bigger than a two pixel distance 

W and the shoulder region is only 1-2 pixels in length. Thus, the direction in which to move the 

CDS boundary is selected for computationally efficiency. 

ffl A second possibility for resetting the pixel values is illustrated by the dashed line 122' in 

Q FIG. 6B. In this case, the values of pixels PI, P2 are set to bowel values. 

E§0 Although the above DBS processes have been described in conjunction with two- 

U dimensional images, it should be appreciated that the same concepts apply equally well to three- 
dimensional images. That is the concepts can be applied along any axis of an image. 

Referring now to FIG. 7, a process for automated polyp detection is described. It should 
2 5 be noted that prior to the processing performed in conjunction with FIG. 7, a polyp template is 
formed. The polyp template is formed by obtaining a CT image which includes a polyp and 
excising the polyp from the image. That excised polyp image is then used as the template. 
Alternatively, one can empirically generate a family of templates that resemble the known 
morphology and density of polyps (or other structure sought to be detected). 
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Processing begins with the step 130 in which an image portion is selected. Processing ten 
proceeds to step 132 in which a polyp identification function is applied to the selected image 
portion. Next, in step 13 4, a correlation is performed between the output of the identification 
function and the image. 

The correlation is performed as: 

[Correlation Matrix] = 2d inverse Fourier Transform of [2d inverse Fourier Transform(T) 
*2d inverse Fourier Transform(I)] 

in which: 

T = the polyp template; and 
I = the image. 

It should be appreciated that the template matching technique can also be performed by 
taking a spherical template and applying the spherical template to voxels (i.e. 3D picture 
elements) formed from a series of two-dimensional images by interpolating values between the 
images as is known. Also, the original template must be transposed to form an operational 
molecule. The operational molecule is the template rotated 180 degrees used as a temporary 
computational transform of the template for the purposes of identifying regions of similarity 
between the template and the image corresponding to the search for polyps within the image. 

It should be noted that if the template is formed from a polyp of a certain size, then the 
template must scaled to detect polyps of different sizes. This can by done by utilizing a four 
dimensional process in which the fourth dimension is a scaling of the polyp template. It should 
also be noted that it may be necessary to modify the polyp template to detect polyps having 
shapes which are different from spheres. 

Referring now to FIG. 8, a second process for automatic detection of polyps, referred to 
as the so-called "rolling ball" technique begins with the step of segmenting the image as shown in 
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step 136. In this step all of the image information is taken and certain portions of the image are 
extracted for use in the analysis. A conventional techniques including but not limited to the so- 
called "marching cubes" technique can be used in the segmentation step. It should be appreciated 
that that the segmentation step and all of the steps show in FIG. 8 can be performed on images 
5 which have been cleansed using the above-described DBSP process. Alternatively, the 

segmentation step and all of the steps show in FIG. 8 can be performed on images provided using 
conventional techniques (e.g. images which have not been cleansed using the DBSP technique). 

Next as show in step 138, a virtual ball is rolled along the bowel boundary. When the ball 
10 is rolled around the perimeter, the system detects changes in direction of the ball to identify 

polyps as shown in step 140. In general, the process of classifying the features of the bowel are 
accomplished by marking three points where turns exist. The points are selected by looking at the 
? changes in slope (e.g. the derivative). Next the distances between the marked points are 
CO computed and the ratios between the distances are used to identify bowel features or 
5j5 characteristics of the bowel. The features of the bowel perimeter are thus classified based upon 
5 the path of the ball and an advance knowledge of the geometric characteristic of the polyp or 
m other irregularity being detected is therefore required. 

£ An example of the rolling ball polyp detection technique is shown in Figs. 8A-8G. 

C| 0 Referring first to Fig. 8 A, an image 142 includes several regions 143a - 143f. As shown in Fig. 
U 8B, region includes a features 150 and 152. A rolling ball polyp detection process is first 
performed on feature 150 as shown in Fig. 8C. 

As shown in Fig, 8C, a test element 154 having the shape of circle (or ball) is moved or 
2 5 rolled along a surface 1 56 which may correspond for example to the surface of a bowel wall. In 
this example, the ball moves in the lumen region 155. It should be appreciated, however, that in 
alternate embodiments, the test element 154 may be provided having a shape other than a circular 
shape. It should also be appreciated that in some embodiments it may be desirable to move the 
ball 154 in a region other than the lumen region 155. When the ball 154 reaches a section of the 
30 surface 156 at a location where two portions of the ball 154 contact two portions of the of the 



surface 156 then a first point 158 is marked. Point 158a corresponds to the point where the slope 
of the surface 156 changes sign. 

After point 158a is marked then the ball continues its path and points 158b, 158c are 
marked. Thus three points are marked where turns exist. Each of the points 158a - 158c are 
selected by looking at the changes in slope (e.g. the derivative). Next the distances between the 
marked points are computed. Here the distance between points 158a and 158b is designated as 
159a, the distance between points 158b and 158c is designated as 159b and the distance between 
points 158a and 158c is designated as 159c. The ratios between the distances 159a - 159c are 
then used to identify bowel features or characteristics of the bowel. In Fig. 8C, the ration is 
computed as follows: the sum of the lengths 1 59a and 1 59b divided by the length 1 53c. The 
ratio of the lengths 159a - 159c along the path defined by points 158a - 158c indicate that the 
structure 152 has the shape of a polyp. This can be accomplished by comparing the ratio value to 
a predetermined threshold value. Thus, the above techniques computes the ratio formed by the 
deflection and travel paths of the test element and uses this information (e.g. by comparing the 
ratio value to a threshold value) to characterize the bowel structure. 

As shown in Fig. 8D, the point 158 is defined by the intersection of lines 159a, 159b. 
Each of lines 159a, 159b project from the center of the ball 154 to the point at which the 
respective surfaces of the ball contact the two points of the surface 156. 

It should be noted that it is important to distinguish collision points (e.g. points 158a, 
158c) from turnaround points (e.g. point 158b). 

Referring now to Fig. 8E, a semi-circular shaped structure 153 exists on the side of a 
triangular shaped structure 152. It should be appreciated that in this case, there are two ratios to 
compute. One ratio value for the structure 152 (i.e. the structure having the triangular shape) is 
formed by the sum of the lengths 161a - 161e divided by the length 161f. Another ratio value for 
the semi-circular shaped structure 153 is formed by the sum of the lengths 161c - 161c divided by 
the length 161g. 
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Referring now to Figs. 8F - 8H, surface 159 forms a structure 160. The rolling ball 
technique is used to define three points 162a, 162b, 162c which define the structure 160. In Fig, 
8F, an average center point 164 has been computed (e.g. via the segmentation technique) and can 
5 be used to provide sufficient information with respect to "turning in" versus "turning away." This 
is accomplished by evaluating the absolute value of slope from collision to collision slope 
(referred to hereinafter as the pre-slope). The absolute value of the preslope will change as the 
ball test point moves around the bowel circumference. When a collision is found, the preslope 
value is used to assess the turn around point(s) and returns to the bowel wall. Turn around points 
1 0 correspond to maxima and minima of the slope encountered by the ball center point as this point 
moves from wall collision to wall collision. 

5 A collision can be defined as when the perimeter of the test ball contacts two portions of 

01 the wall as illustrated in Fig. 8D. By setting the perimeter of ball correctly - as determined by 
"fy> empirical evaluation - one can select feature changes of a polyp or fold from the background 
'M concavity of the bowel perimeter. 

□ It should be appreciated that the rolling ball technique can also be performed in three 

f dimensions (i.e. both around a circumference and up and down a longitudinal axis of a short tube) 

p 0 as shown in FIG. 81. In Fig. 81, a series of axial slices 1 65a - 1 65c are used to construct a three- 

2 dimensional image of bowel sections 166 (small bowel loop) and 167 (loop of bowel). Regions 
168a - 168e correspond to the portions of the bowel which would be present in the axial images 
(e.g. as shown in Figs. 1 A and 1C). In the case where the rolling ball technique is applied to a 
three-dimensional bowel image, care must be taken to insure that each region of the bowel "tube" 

25 is examined. 

Referring now to Fig. 9, a third process for automatic detection of polyps, referred to as 
the so-called "distance search" technique begins with the step of selecting a polyp template 170. 
The template is selected having a predetermined shape. The shape is preferably selected to 
3 0 correspond to the shape of the bowel anomaly. 



Referring briefly to FIG. 9 A, for example, a template 180 includes a local window 
boundary 182 and a template pattern 184. The template pattern is here selected to be a circle 184 
since the polyps tend to have substantially circular shapes as seen in 2D images. If applied in 3D, 
then spherical, or semi-sphere or other forms could be used by utilizing a 3D implementation. It 
should, however, be appreciated that in those applications in which the lesion sought to be 
detected had other than a circular shape, the template pattern would be selected accordingly. For 
example, if the lesion sought to be detected had a substantially triangular shape, then the template 
pattern would also be selected having a substantially triangular shape. The side 180a of the local 
window 180 is selected having a length which is as large as the largest lesion to be detected. The 
template pattern 184 is selected having a dimension as small as the smallest lesion to be detected. 

Referring again to FIG, 9, once the window size and template pattern and size are 
selected, the local window 182 is scanned across the image of the bowel perimeter once it has 
already been segmented. Next, as shown in step 174 in FIG 9, the distances between points on 
the template 184 and the perimeter points within the window 182 are made. 

This may be more clearly understood with reference to FIGs. 9B - 9D, in which an image 
186 is segmented to include only the bowel perimeter 188. The bowel perimeter includes a fold 
190 and a lesion 192. The window 182 is placed over the image 186 and moved across the entire 
image 186. When the window 182 reaches a location in the image 186 in which a portion of the 
bowel boundary is within the window and within the template 184 (as shown in FIGs. 9C, 9D) the 
distance between the points on the template 184 and the perimeter points within the window 182 
are made. In one embodiment, the distances between 10-30 test points on the unit test circle and 
boundary points present within the moving frame should initially be made. Those of ordinary skill 
in the art should appreciate of course that the precise number of points used is not critical and that 
in some applications, it may be desirable or necessary to use a number of points fewer than 10 or 
greater than 30. 



Referring now to step 176 and FIG. 9E, a determination is made as to whether the 
distances are equal. When the test circle becomes centered within a circular lesion, then the 
distances between the circle perimeter and boundary points of the lesion becomes equal and the 
standard deviation of a group of those distances approaches a minimum value (See Figs. 9F and 
9G). Its should be noted that the above process is carried out in one plane and it is necessary to 
search in three orthogonal planes to be complete. The point where the standard deviation of a 
group of those distances approaches a minimum value is the point which should be marked as a 
center of a suspected lesion (i.e. identified as a region of interest as indicated in step 178). After 
the technique is run in all three planes, those lesions that were tagged in 2 of 3 or 3 of 3 planes 
can be finally tagged as suspicious regions. This technique can thus be used to distinguish a fold 
from a polyp. It should be noted that depending upon the feature to be selected, other statistical 
tools, aside from standard deviation can be used to identify regions of match between the roving 
template and features of the segmented image. 

The technique described above in conjunction with FIGs. 9-9G,thus searches for patterns 
of distance. The technique is thus relatively computationally expensive but it is also relatively 
rigorous. One advantage of this technique is that t does not matter how large the template is 
relative to the polyp so one scanning works for all lesions. 

In Fig. 9G, a plot 209 of distance around the template vs. template position is shown. 
Region 210 of curve 209 is relatively flat thus indicating that a region of the image contained a 
shape which matched the shape of the template. When the template is provided having a round 
shape (e.g. as shown in Figs. 9E, 9F), the flat region 210 indicates that a corresponding shape is 
found in the structure being searched. 

All references cited herein are hereby incorporated herein by reference in their entirety. 

Having described preferred embodiments of the invention, it will now become apparent to 
one of ordinary skill in the art that other embodiments incorporating their concepts may be used. 
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It is felt therefore that these embodiments should not be limited to disclosed embodiments, but 
rather should be limited only by the spirit and scope of the appended claims. 

What is claimed is: 

5 
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